DOE/ER/40427-28-N96 



On 
On 



C/3 



The role of chiral symmetry in two-pion exchange nuclear potential* 

C.A. da Rochal 

University of Washington, Dept. of Physics, Box 351560, Seattle, Washington, 98195-1560 

M.R. Robilotta* 

University of Washington, Dept. of Physics, Box 351560, Seattle, Washington, 98195-1560 

and 

Instituto de Fisica, Universidade de Sao Paulo, CP. 20516, 01452-990, Sao Paulo, SP, Brazil 

J.L. Ballot 

Division de Physique Theorique, Institut de Physique Nucleaire, F91406 Orsay Cedex, France 

(August 96) 



We evaluate the two pion exchange contribution to the nucleon-nucleon potential in configuration space using a minimal chiral 
model containing only pions and nucleons. We argue that this model has nowadays a rather firm conceptual basis, which 
entitles it to become a standard ingredient of any modern potential. The main features of this model is that the scalar-isoscalar 
Oh- component of the interaction is relatively small, due to cancellations between large terms, and fails to reproduce the intermediate 
range attraction in the central channel. We show that chiral symmetry is the responsible for these large cancellations in the 
two-pion exchange nucleon-nucleon interaction, which are similar to those occuring in free pion-nucleon scattering. Another 
' feature of the model is that these results do not depend on how chiral symmetry is implemented. 
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I. INTRODUCTION 

The long range part of the nucleon-nucleon (NN) in- 
teraction is ascribed to the exchange of one pion (OPEP) 
and is very well known. The medium range region, on the 
other hand, is much more controversial, since in the liter- 
ature one finds various competing theoretical approaches, 
as dispersion relations, field theory or just phenomenol- 
ogy. In all cases, the medium range interaction is asso- 
ciated with the exchange of two pions (TPEP), whereas 
there is a wide variation in the way short distance ef- 
fects are treated. In order to reproduce the data, pa- 
rameters are used which either reflect knowledge about 
other physical processes or are adjusted ad hoc. The 
differences in all these approaches are enhanced in the 
isospin-symmetric nuclear matter, where the OPEP van- 
ishes and one must know the asymptotic behavior of the 
TPEP. 

The role of the TPEP in the framework of chiral sym- 
metry has recently attracted considerable attention, es- 
pecially as far as the restricted pion-nucleon sector is con- 
cerned [E]@J3M|tfiJ. It is well known that this process is 



* Complete version of the contribution presented at the 14 th 
International Conference on Particles and Nuclei (PANIC 96), 
Williamsburg, VA, 22-28 May 1996. 

'Fellow from CNPq, Brazilian Agency. Eletronic address: 
carocha@phys.washington.edu 

■''Eletronic address: robilotta@if.usp.br 

^ Unite de Recherche des Universites Paris 6 et Paris 11 as- 
sociee au CNRS. Eletronic address: ballot@frcpnll.in2p3.fr 



closely related to the pion-nucleon scattering amplitude 



Chiral symmetry at hadron level may be implemented 
by means of either linear or non-linear Lagrangians. The 
fact that no serious candidate has been found for the a 
meson favours the latter type of approach. There are 
two forms of non-linear Lagrangians which are especially 
suited for the 7rN system. One of them is based on a 
pscudoscalar (PS) 7rN coupling supplemented by a scalar 
(S) interaction, equivalent to the exchange of an infinitely 
massive a meson, and denoted as PS + S scheme (Eq. [l]) . 
The other one employs a pseudovector (PV) 7rN interac- 
tion and a vector (V) term, which could represent the 
exchange of an infinitely massive p meson, constituting 
the PV + V scheme (Eq. §). 



Cps+s = gip [V fl - 4> 2 + it ■ 075 ip , (1) 
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Both approaches yield the very same amplitude for the 
7rN scattering for the isospin-symmetric amplitude. The 
fact that physical results should be independent of the 
representation used to implement chiral symmetry was 
discussed in very general terms by Coleman, Wess, and 
Zumino §. 

As far as nucleon-nucleon scattering is concerned, the 
two-pion exchange amplitude up to order 1//^ is given 
by five diagrams, usually named box (□), crossed box 
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FIG. 1. Loop diagrams for the two pion exchange NN potential calculated in the minimal chiral model. 



(X), triangle (A + V) and bubble (0), given in Fig. |l|, which constitute the minimal chiral model or the TPEP. 
The first two diagrams contain only nucleon propagators and are independent of chiral symmetry, whereas the triangles 
and the bubble involve the scalar interaction and hence are due to the symmetry. When one considers the potential 
instead of the amplitude, the iterated OPEP has to be subtracted from the box diagram. 



II. PARAMETRIZATION OF THE TPEP 



Our calculation of the TPEP is based on the Blankcnbecler-Sugar reduction of the Bethe-Salpeter equation. Its 
dynamical content is associated with the five diagrams displayed in Fig. |l| Therefore we label the corresponding 
individual contributions by (j , A , X and □, where the last one also includes the subtraction of the iterated OPEP. It 
has the general form 
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where the spin operators are given by Oss — -cr^ 2 \ Ols 
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whereas erW and represent spin and isospin matrices for nucleon (i). 

In the case of the bubble diagram, the leading contribution to the asymptotic potential can be calculated analyti- 
cally ||; its result sets the pattern for the parametrization of the other components of the force. 

Our numerical expressions represent the various components of the potential in MeV, and are given in terms of 
the adimensional variable x = /ir, where [i is the pion mass. We keep the axial 7rN coupling constant gA as a free 
parameter and adopt the values /i = 137.29 MeV and f„ = 93 MeV for the pion mass and the pion decay constant 
respectively. 

In general, the parametrized expressions reproduce quite well the numerical results for the complete theoretical 
calculations, given in Ref. Q, except for a few cases and regions where the discrepancies become of the order of 
0.25%. Our results are listed below. 



A. Central Potential 



The profile function for the central potential has the following common multiplicative expression 

V 2U J X 2 y/X 



The parametrization of each diagram gives 

3 6.54068 1.26190 0.130706 
— i 1 
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V ( f(x) = F c (x) -275.364 — + — — + —— (5) 
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B. Spin-spin potential 

The multiplicative factor for the spin-spin potential is the same as that of the central potential, and receives 
contributions from the box and crossed diagrams only, which are given by 

VSS( X ) = F c (x) (0.408084+ + - - 0.215829 • e— } (9) 

V*"(x) = F c (x) (0.399845 + + °^02 - °-^± + 0.0371333 * • e— j (1Q) 



C. Spin-orbit Potential 

The spin-orbit multiplicative function is 



and individual contributions are: 



TrL3 , . ^ , . f r 00 _ 5.51078 0.994157 0.217562 0.0336067 0.002446201 

V LS (x) = F LS (x) |-5.88744 — - + —g -j— + — -j -j f (12) 
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V^z) = F LS (x) |7.34548 + - - 3 ^ 025 + (7.48798 - 0.448484 • x + 0.391431 • x 2 ) ■ e -"-4i9984.* j (13) 

T *, . V£ s (a;) f 5.67235 10.6417 14.7389 12.6506 6.27374 1.66637 0.1842641 , x 



(15) 
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V^ s {x) = F LS (x) I -1.19527 + - - \ ■ [1 - 0.0160259 • x] 1 

I X X X J 

D. Tensor Potential 

The common factor for the tensor potential is 

It receives contributions from the box and crossed diagrams only, which have the form 
V£( X ) = F T (x) {-0.204041 - Q" 51 " 720 + + (1.32932 - 0.939553 • x + 0.706050 • x 2 ) ■ e -2.29686.,| 

TrT . . ^ , ,[ 0.521123 0.352463 0.135028 0.0303012 0.002978401 

V T(x) = F T (x) |-0.246349 — + — ^ — + — ^ ^ j (18) 
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FIG. 2. Profile functions for the bubble (Q) and triangle 
(V) scalar- isoscalar potentials and for their sum (S), showing 
a strong cancellation between these two contributions. The 
graph at right is an amplification. 



FIG. 3. Profile functions for the chiral (S) and pseudoscalar 
(PS) scalar- isoscalar potentials and for their sum (M), show- 
ing another strong cancellation between these two contribu- 
tions. The graph at right is an amplification. 



III. RESULTS 

Using the minimal chiral model potential described in 
last section we evaluated the scalar-isoscalar component 
of the TPEP and the phase shifts for some singlet waves. 
We adopted gA = 1-33 @. It is possible to notice two 
important cancellations within the scalar-isoscalar sector 
of the 7T7rE-NNP. The first of them happens between the 
triangle and bubble contribution, as shown in Fig. |2[ 

The other one occurs when the remainder from the 
previous cancellation (S) is added to the sum of the box 
and crossed box diagrams (PS). In this last case, the di- 
rect inspection of the profile functions for the potential, 
given in Fig |^, provides just a rough estimate of the im- 
portance of the cancellation, since the iterated OPEP is 
not included there. 

The second cancellation can be better studied in the 
NN scattering problem, since the amplitudes obtained 
by solving a dynamical equation include automatically 
the iterated OPEP and depend very little on the way 
the potential is defined. In our case, we are interested in 
exhibiting the effects associated with chiral symmetry in 
the scalar-isoscalar two-pion exchange channel. There- 
fore we concentrate our study on singlet channels, where 
the strong effects associated with tensor OPEP interac- 
tions are not present. A suitable choice of observables 
also allow a separation of the dynamical effects according 
to their range. This is particularly useful in this problem, 
since it is well known that the use of the Chiral Perturba- 
tion Theory (CHPT) is associated with the inclusion of 
undetermined counterterms in the Lagrangian, involving 



higher orders of the relevant momenta (l]|||ic]]. How- 
ever, in configuration space, these counterterms become 
delta functions which affect just the origin and hence are 
effective only for waves with low orbital angular momen- 
tum. In our derivation of the NN potential we used a La- 
grangian which did not contain these contact terms and 
hence it is suited for medium and long distances. Thus, 
in order to avoid these undetermined short range effects, 
we consider only the 1 Z?2, G4, Fa, and 1 H$ waves. For 
each channel, we decompose the full NN potential V as 

V = U n + U PS + U s + U c , (19) 

where U-n is the OPEP, Uc represents the short ranged 
core contributions, Ups is due to the box and crossed box 
diagrams whereas Us is associated with the chiral trian- 
gle and bubble interactions. Using the variable phase 
method, it is possible to write the phase shift for angular 
momentum I as |ll| , ^2| 

/>oo 

h = ~\ drVPj. (20) 
k Jo 

In this expression, the structure function Pi is given by 

Pi = ]e cos Di — hg. smDi , (21) 

where ji and hi are the usual Bessel and Neumann func- 
tions multiplied by their arguments and Di is the variable 
phase. Using the decomposition of the potential given in 
Eq. ||, one writes the perturbative result 
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FIG. 4. Contributions for the long-OPEP (ttL), iterated 
OPEP (ttI), pseudoscalar (PS) and chiral (S) terms of the 
potential to the phase shifts for 1 Z>2 and 1 Gi waves. The total 
phase shifts are indicated by (T). The experimental result is 
given by the dashed line. 
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FIG. 5. Contributions for the long-OPEP (wL), iterated 
OPEP (ttI), pseudoscalar (PS) and chiral (S) terms of the 
potential to the phase shifts for 1 Fs and 1 H$ waves. The total 
phase shifts are indicated by (T). The experimental result is 
given by the dashed line. 
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In this expression, the first term represents the pertur- 
bative long range OPEP (ttL), the second the iterated 
OPEP (nl) , the third the part due to the box and crossed- 
box diagrams (PS), the fourth the contribution from chi- 
ral symmetry (S). The last one is due to the core and 
vanishes for waves with I ^ 0. 

In Fig. [| we show the partial contributions of the sin- 
glet even waves 1 I?2 and 1 Gi to the phase shifts as func- 
tions of energy. 

Observing Fig. ||, one sees that the minimal chiral model 
with a core added fails to reproduce the energy depen- 
dence of 1 £>2 wave and is resonable for the 1 Gn wave. 
This not-so-well agreement was already expected since 
the minimal chiral model does not include the complete 
dynamics of the ttN scattering. Mesons and baryons res- 
onances like p and A are very important for higher ener- 
gies and lower orbital angular momenta £. As one goes 
to higher values of £, the role of the one and two pion 
exchange becomes more important. This can be seen in 
Fig. |5[ where we show the partial contributions of the 
singlet odd waves ^-F^ and 1 H§ to the phase shifts as 
functions of energy. 

A remarkable feature of these two waves is that the net 
result is given just by the long-OPEP contribution, mean- 
ing that the medium range contributions cancel entirely. 



In these channels, the iterated OPEP is noticeable and 
the relationship 5t)„i + Se)ps = —Si)s holds. This is an 
important feature of the chiral symmetry and explains 
why the first theoretical models for the TPEP in the 50 's, 
which do not have the S term, spoiled the good agree- 
ment achieved by the OPEP alone for phase shifts with 
large I. Moreover, it is possible to see another impor- 
tant features of the TPEP, namely the iterated OPEP 
contribution is comparatively small, indicating that am- 
biguities in the definition of the potential do not have 
numerical significance. 

Our results show that chiral symmetry, in the re- 
stricted pion-nucleon sector, is responsible for large can- 
cellations in the two-pion exchange interaction. This pro- 
cess is therefore similar to threshold pion-nucleon or pion- 
deuteron |l3| scattering amplitudes, where the main role 
of the symmetry is to set the scale to the problem to be 
small. 



IV. PERSPECTIVES 

A shortcoming of the minimal chiral model is that it 
fails to reproduce experimental information in the case of 
the intermediate irN amplitude that enters the TPEP. In 
order to overcome this difficulty, one may extend the ap- 
proach so as to encompass other degrees of freedom. This 
possibility was recently considered by the introduction 
of the A resonance , which improved considerably the 
predictive power of the method. Another way to achieve 
this goal is to introduce the empirical information that 
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FIG. 6. Diagrams associated with the ttN amplitude, (a), (b), and (c) represent the minimal chiral model, whereas (d) 
corresponds to the net effect of the HJS coeficients and is denoted by R, for "rest" . 



is missing in the intermediate ttN amplitude in a model independent way, with the help of the Holder, Jacob and 



Strauss |14 T^Jl^l subthreshold coeficient, as proposed by Robilotta M. The detailed knowledge of the ttN amplitude 
provided by these coeficients allows a precise determination of the missing part of the TPEP. as illustrated in Fig. ||. 
This approach is very similar to in spirit to that followed long ago by Tarrach and Ericson jL7j , who explored in detail 
the analogy between the TPEP and Van der Waals force. 

A detailed calculation for the leading term in this approach was already done |^| , and a complete calculation will be 
presented soon lis). 



ACKNOWLEDGMENTS 



We thank Prof. U. van Kolck for helpful conversations about Chiral Perturbation Theory and the NN force. This 
work was partially supported by U.S. Department of Energy. The work of C.A. da Rocha was supported by CNPq 
Brazilian Agency. 



[1] C. Ordonez and U. Van Kolck, Phys. Lett. B 291, 459 (1992). 

[2] L.S. Celenza, A. Pantziris and CM. Shakin, Phys. Rev. C 46, 2213 (1992). 

[3] J.L. Friar and S.A. Coon, Phys. Rev. C 49, 1272 (1994). 

[4] C.A. da Rocha and M.R. Robilotta, Phys. Rev. C 49, 1818 (1994). 

[5] M.C. Birse, Phys. Rev. C 49, 2212 (1994). 

[6] C. Ordonez, L. Ray, and U. Van Kolck, Phys. Rev. Lett. 72, 1982 (1994). 
[7] C.A. da Rocha and M.R. Robilotta, Phys. Rev. C 52, 531 (1995) 

[8] S. Coleman, J. Wess, and B. Zumino, Phys. Rev. 177, 2239 (1969); C.G. Callan, S. Coleman, J. Wess, and B. Zumino, Phys. 

Rev. 177, 2247 (1969). 
[9] M.R. Robilotta, Nucl. Phys. A 595, 171 (1995). 
[10] V. Bernard, N. Kaiser, and Ulf-G. Meissner, Int. J. Mod. Phys. E 4, 193 (1995). 

[11] F. Calogero, Nuovo Cimento 27, 261 (1963); Variable Phase Approach to Potential Scattering, (Academic Press, New York, 
1967). 

[12] J.L. Ballot and M.R. Robilotta, J. Phys. G 20, 1595 (1994). 

[13] M.R. Robilotta and C. Wilkin, J. Phys. G 4, L115 (1978). 

[14] G. Hohler, H.P. Jacob, and R. Strauss, Nucl. Phys. B 39, 237 (1972). 

[15] R. Kock and E. Pietarinen, Nucl. Phys. A 336, 331 (1980). 

[16] G. Hohler, Landolt-Bornstein Numerical Data and Functional Relationships in Science and Technology, group I, Vol.9, 

subvol. B, part 2, Ed. H. Schopper, (1983). 
[17] R. Tarrach and M. Ericson, Nucl. Phys. A 294, 417 (1978). 
[18] M.R. Robilotta and C.A. da Rocha, to be published. 



G 



